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Abstract. We search for Borromean three-body systems of identical bosons in two dimensional geometry, 
i.e. we search for bound three-boson system without bound two-body subsystems. Unlike three spatial 
dimensions, in two-dimensional geometry the two- and three-body thresholds often coincide ruling out 
Borromean systems. We show that Borromean states can only appear for potentials with substantial 
attractive and repulsive parts. Borromean states are most easily found when a barrier is present outside 
an attractive pocket. Extensive numerical search did not reveal Borromean states for potentials without 
an outside barrier. We outline possible experimental setups to observe Borromean systems in two spatial 
dimensions. 
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1 Introduction 

In quantum mechanics two particles in three spatial di- 
mensions (3D) can attract each other without forming a 
bound state. However adding a third particle can make a 
three-body system bound. Such bound three-body struc- 
tures, where each pair of particles is unbound, are called 
Borromean states. They turn out to be rather abundant 
in subatomic physics where many examples are found and 
studied, see e.g. [JJEJ- Those examples raised the question 
about properties of the potentials that determine the pos- 
sibility for Borromean binding. One can trace this discus- 
sion for three-dimensional geometry in numerous papers, 
e.g. [210 • It turns out that in 3D finite range potentials of 
the form gV{r) most likely have a region of the parameter 
g where three, but not two particles are bound. 

In two spatial dimensions (2D), however, this question 
is not well-established. The result has to be different, be- 
cause two-body binding is easier to achieve in 2D than in 
3D. This is shown already for two particles with gV < 0. 
Such potentials support two- and three-body bound states 
even when g — > [5] . The thresholds for binding of two and 
three-body systems are identical and Borromean struc- 
tures cannot exist. Moreover, numerical investigations in 
momentum space [BUZ] and coordinate space [Sl CTTUlfTTl 
n"2"l[T5] revealed that in 2D some other classes of poten- 
tials have no Borromean states. 

Previous investigations also show, that it is not entirely 
impossible to have Borromean systems in 2D. However, to 
the best of our knowledge, only one example of an appro- 
priate two-body potential can be found in the literature 
[TO] . The purpose of this report is to identify conditions 
for the two-body potentials that can produce Borromean 
systems of three identical bosons in two dimensions. We 



believe that using advanced experimental techniques it is 
possible to create setups where these conditions are satis- 
fied. This introduces structures with new few- and many- 
body properties. This investigation also puts limitations 
for a zero-range formalism that is used widely to describe 
few-body dynamics in 2D [T4lll5llT6] . This formalism is a 
powerful tool which has been proven to give correct re- 
sults for weakly-bound systems in 2D that have identi- 
cal three and two body thresholds. However it can not 
describe Borromean states where finite-range techniques 
have to be used. 

We believe that this investigation is timely and impor- 
tant. First, because the experimental techniques in cold 
atomic gases are rapidly improving in 3D and are cur- 
rently being adjusted to 2D. Second, Borromean system 
has three particles which strongly suggests that traditional 
many-body correlations built on two-body properties now 
may turn out to be completely different if three particles 
are considered as building blocks. Our results might pro- 
vide guidence to the tuning of the potentials in order to 
achieve new many-body structures. 

In this report we briefly discuss conditions for two- 
body binding in section II. We search for characteristic 
properties, e.g. shape of the potential V and strength g, 
which allow binding of three but not of two particles in 
section III. We ouline experimental setups to observe Bor- 
romean states in section IV. Finally, in section V we briefly 
summarize and conclude. 



2 Two-body problem in two dimensions 

According to the definition Borromean states, if exist, ap- 
pear for two-body potentials at the edge of two-body bind- 
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ing. An additional particle then provides the glue to form a 
three-body bound state. We therefore first must establish 
the criteria for a two-body potential to support a bound 
state [T7J. This is known in the limit of weak potentials, 
g — > 0, where the net volume, g J Vrdr, simply has to 
be non-positive. However, the existence of bound states 
for stronger potentials is not solely determined by the net 
volume. Thus we reformulate the two-body problem with 
the aim of extracting suitable criteria for binding. The 
qualitative behavior is supplemented by a full quantita- 
tive analysis for simple potentials. 



2.1 General properties 

The radial Schrodinger equation for two particles in 2D is 

(E 2 - 9 V)4>m , (1) 



— (--— r— + —\d 

2fi\ r dr dr r 2 



where r is the relative coordinate, E 2 the two-body en- 
ergy, gV(r) the potential, \i the reduced mass, and AI = 
0,1,2,... is the angular quantum number, <pM the wave 
function of the relative motion. For spherical potentials 
M is conserved and M = therefore characterizes the 
ground state for identical bosons. We shall in the follow- 
ing only consider M = and omit any related index. 

We demand the wave function, <j)(k, r), to be regular in 
zero, through <f>(k, r = 0) = 1, where k is the wave number 
given by E 2 = H 2 k 2 /(2(i). The condition for a bound state 
can then be written as [P8l[T9l 



■ poo 

1 + g— / drr<j>(k,r)V(r)H^ } (kr) = , (2) 
2 Jo 

where Hq is the first Hankel function of order zero. The 
solutions correspond to bound states with k = i\k\. We 
focus on the weak binding limit where k approaches zero 
and eq. © reduces to 



i 



g / drr0(O, r)V{r)\i- 



where 7 is Euler's constant and the regular zero energy 
solution to eq. (HJ is 



(4) 



0(0, r) = 1 + g / dss In (r/s) V(s)0(O, s) . 
Jo 

Then eq. can only be fulfilled for small k when 

/>oo 

g / drr<j){0,r)V{r) = . 



(5) 



Eqs. and (|4]) define appearence of the bound state 
in 2D. This result might also be obtained using the Jost 
function formalism, see [2"0] . 

From eq. Jl} we get immediately 



" R h 2 ( d 

drr<f>{0,r)V{r) = _ (r— 0(0, r) 



0, (6) 



where the vanishing result is obtained for finite-range po- 
tentials where V(r > R) = 0. Thus the slope of the zero- 
energy wave function is zero at r = R. 

It is obvious that for purely attractive potentials, gV < 
0, the condition from cq. ([S]) can be achieved only for 
g = 0. This is a consequence of the fact that purely at- 
tractive potentials always provide at least one two-body 
bound state in 2D [S] . In this case the thresholds for bind- 
ing of two and three-body systems has to coincide hence 
Borromean states are ruled out [10 . 

It gives us the first necessary condition for Borromean 
binding: the potential, gV, has to contain both positive 
and negative parts. We shall proceed by defining shapes 
V± and dimcnsionless strengths A± > of the positive 
and negative parts. The total potential is then given by 
gV(r) — A + y + + A_y_. Binding is subsequently achieved 
by sufficient increase of A_, and binding is reduced by 
increase of A+. We define the critical repulsive interaction, 
A^ r) (A_), through eq. © by 



drr(j>{Q 1 r){\ c lV+ + X-V-) = 



(7) 



This implies that a Borromean system can only appear 
when A+ is larger than AST, since two particles are un- 
bound in this regime. Hence a Borromean system for given 
A_ might be found in an interval where A+ is larger than 
AST, although perhaps an interval of very limited exten- 
sion. The solution of eq. (@| and the definition in cq. © 
then provide crucial information about the most likely re- 
gion for occurrence of Borromean systems. 

Let us use eq. ([7]) for potentials, where knowledge of 
the wave function is not needed. First we show the well 
established result that any weak potential with negative 
or zero net volume has at least one bound state (51 IT8l[2Tj . 
When the attractive part is very weak we can appoximate 
the wave function in eq. ((U) by the first term. This gives 



A^ r 



-A- 



J drrV- (r) 
J drrV+(r) 



+ 0(e --). 



(8) 



j (3) Then from eq. ((Tj) we see that when A+ < \+ the net 



r=R 



volume is less than 0(e A - ) <C A_ and at least one 
bound state is present, whereas the system is unbound 
for A+ > AST where the net volume is positive. This result 
leads us to a second necessary condition for a Borromean 
system, that is g J Vrdr > 0, as only those potentials have 
a region of g without two-body states. Second, we consider 
the two delta-shell potentials, gV = K 2 j (2{id 2 )(A + 5(r jc — 
1) — XS(r/d— 1)), which is infinitely large for a given ra- 
dial distance and zero otherwise. In the limit when c — > d 
we immediately conclude from cq. ([JJ and continuity of 
the wave function that A^T = — A_. 

Other potentials localize the wave function in the at- 
tractive region, which increases J drrV-(r)<f>(Q,r) and de- 
creases / drrV + (r)<p(0, r). Eq. © then strongly suggests 

dX cr 

a steep increase of at a given sufficiently large value 
of A_. Or, in another words, if the potential is of finite 
range then for A^jT — > 00 a finite Al r exists that will bind 
two-body system. 
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2.2 Square well with barrier or core 

The threshold conditions are not easily derived for ar- 
bitrary potentials with sizable attraction and repulsion. 
However, the general equations, eqs. © and ([1]), are di- 
rectly applicable for simple potentials where zero-energy 
wave functions are known. From properly selected poten- 
tials we can then extract correct qualitative properties for 
more general potentials. We choose to study a solvable 
model containing all the crucial features, that is 



V = X-V- + X+V+ 



2fiR 2 s 

x + h 2 

2/j.Rf 





r < R s 

R s <r<Ri 

r>Ri, 



(9) 



where R s and Ri are positive radii. This potential is useful 
for describing qualitative features of more general poten- 
tials, see e.g. [22j. The zero-energy solution then has the 
form 

AJaihr) r<R s 
0(0, r) = { BI (k 2 r) + CK (k 2 r) R s < r < Ri (10) 
DK (k 3 r) r>Ri, 

where fci = \fX_/R s , k 2 = \fX+/R s , k 3 0, J ,Io,K 
are Bessel functions with the usual definitions [25] . and 
the constants B, C, D are defined through A by matching 
at the points r = R s and i?;. One of these equations is 
the quantization condition providing the energy E. For 
E — we use instead eq. ([7]), which by use of xZq(x) = 
j^xZi(x),Z l = I^Y^Ji, xK {x) = --^xKi(x) and the 
wave function in eq. (|10p can be integrated to give 

ah(k 2 Ri) = bK 1 (k 2 Ri), 
where the constants a and b arc defined by 



(11) 



a = k 2 J (k 1 R s )K 1 {k 2 R s ) - k 1 J 1 (k 1 R s )K {k 2 R s ), 
b = k 2 Jo(k 1 R s )I 1 (k 2 R s ) + kiJi(kiR a )I {k 2 R.). (12) 

The derivation of eqs. ([TT]) and ([T2"]) employed the bound 
state boundary condition and the relations Jq(x)Y\{x) — 
Ji(x)Y {x) = - A and I (x)Ki(x) + Ii(x)K (x) = ±. 

First we start with overall weak potentials, that is 
k 2 R s — > and k\R s — > 0. Then the properties of the 
Bessel functions give a — > l/R s and b — > R s (kf + k\)/2. 
Using that J\(x) ~ x/2, K\(x) ~ 1/x we find that eq. 
is equivalent to / Vrdr = 0, and coinciding with eq. ([8]). 

The other limit of a sizable barrier where k 2 R s ^> 1 
can also be found analytically. After some manipulations 
using properties of the Bessel functions, we conclude that 
Eq. ([TT]) can be approximated by 



hRsJoihRs) - hRshfaRs) = , 



(13) 



where we assumed that Ri ^ R s . Thus, the solution has to 
be near the nodes of the Bessel function, since Jo(kiR s ) = 
k\R s Ji(kiR s )/(k 2 R s ) — > 0. Continuity of the wave func- 
tion at r = R s combined with the exponential decrease 
by penetrating into a barrier then in turn implies that the 
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Fig. 1. The critical strength, AJjT, on the positive part of the 
potential as function of the strength, A_, on the negative part. 
We use R s /Ri = 1/2. The two curves are for the square well 
with a barrier (solid green), eq. (|lip . and the inverted with a 
core and no outer barrier (dashed blue), eq. (|15[) . The vertical 
dot-dashed lines are the corresponding asymptotic values, A!T. 



wave function must be severely diminishing with increas- 
ing barrier. 

The critical value, A^T, are shown in Fig. [JJas function 
of A_, see eq. ([TT]) . The behavior is typical and thus char- 
acterizes also the behavior for less schematic potentials. 
The linear behavior for weak potentials (small A_) is as 
predicted in eq. ([8]). As A_ increases the critical repulsion 
rises steeply and at some point, Al r , it cannot compen- 
sate to avoid a bound state. Larger attraction, A_ > A^ r , 
always provides two-body binding. 

The potential in eq. (O could as well be understood 
with positive and negative parts interchanged. 



V = X-V-+X+V+ 



A_/i 2 



r<R a 
R s <r<Ri 



r>Ri 



(14) 



Then the corresponding condition in eq. (jllj) has to be 
modified as achieved most easily by analytic continuation. 
This implies the use of k\ = R s ,k 2 = ^yxZ/Rg, 

resulting in 



a.J 1 (k 2 Ri) = bY 1 (k 2 R l 



(15) 



where the constants correspondingly have to be changed 
to 

~a = fc 2 / (fc 1 i? s )r 1 (fc 2 i? s ) + hhikMYofaRs), 
b = k 2 I (kxR s ) Ji{k 2 R s ) + kxh(kiR s )Jo{k 2 Rs). (16) 

For weak potentials we again find a —> 1/R S , b —> (k 2 R s + 
k\R s )/2, which as before leads to the condition for binding 
with zero net volume. 

However, these potentials of finite size with k\R s 3> 1 
behave qualitatively different from the potential in eq. ^ . 
In this limit it is again possible to reduce eq. ([T5]) to give 

Jo(*aiOT(*2-Ri) - Ji{k 2 Ri)Y {k 2 R s ) = 0. (17) 
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Fig. 2. The dimensionless volume, h, of the attractive part 
as function of the ratio, s = R s /Ri, of inner and outer radius 
for the square well potential with a repulsive inner core, see 
eq. (JXTJ) . The horizontal dashed line is the volume needed to 
bind two particles when an outer barrier is present, see eq. (I13I) . 



This equation provides the boundary condition for a wave 
function which is zero at r = R s and has zero derivative at 
r = Ri. This condition for binding by an attractive square 
well then depends on both k2R s and A^i?;. These two pa- 
rameters can be expressed in terms of the dimensionless 
measure of the attractive volume, i.e. h = k^(Rf — R 2 ), 
and the ratio between the two radii s = k2R s / {k^Ri) < 1, 
that is k%Rf = h/(l - s 2 ) and k 2 R 2 = hs 2 /(l - s 2 ). Ex- 
pressed by h and s, eq. (fTT|) determines h as function of s 
as shown in Fig. [3] For s — > 1, h ~ and for small val- 
ues of s: h ~ — 1/ In s. We note continuous approach to the 
s = limit which corresponds to the binding in an over- 
all attractive potential where only infinitesimal strength 
or volume is necessary. The other limit of s = 1 exhibits 
that the potential must be infinitely deep when the two 
radii approach each other at a finite value. This variation 
is completely different from the condition of binding for 
the potential from eq. © where the volume for deep po- 
tentials is independent of the barrier dimension. We also 
solve eq. (|15p for given ratio Ri/R s and exhibit A^ as 
function of A_ in Fig. [TJ Again we find slow convergence 
towards the asymptotic value obtained from eq. (fTT)) . 

In ID we consider potentials given by eqs. © and (|14| 
for x > and infinite wall for x < 0. We write then, in- 
stead of eq. ([T3"|) . k\R s = 7r, where only k\R s is present 
and eq. fTT]) is replaced by k2(Ri — R s ) = tt/2, where only 
the difference Ri — R s appears due to translational in- 
variance [5]. In 3D we get a radial equation which is the 
same as in ID. Moreover we also have the same boundary 
conditions, because we assume that ID potentials are infi- 
nite wall for negative x. Thus, in 3D we get the threshold 
conditions, k\R s = 7r and k2{Ri — R s ) = t/2 for potentials 
from eqs. ((9]) and (f!4|) respectively. 

So far we can conclude that only potentials with pos- 
itive volume integral can provide Borromean states. Also 
we see that the strength of the attractive part, A_, can 
not be larger than AfT, as A_ > A1 r always produce a 



two-body bound state. This is the knowledge that we can 
extract without solving three-body problem. 



3 Three particles in 2D 

We consider the three-body Schrodinger equation for three 
identical bosons with A+ larger than A"'. If this system has 
a bound three-body state we call it Borromean. To make 
this procedure formal we define (A_) as the threshold 
value of A + for binding of the three-body system (for a 
given A_ , potentials with A + > can not provide three- 
body bound states). Since > AST, Borromean states 
then appear if A c + ^ A+, for repulsive strengths in the 
interval A"' < A + < A+. We want to establish which 
potentials can provide Borromean binding. 

3.1 Simple cases 

It is entirely possible that Borromean systems do not exist, 
that is A+ = A"' . We shall here provide a number of such 
examples, some of which we have already mentioned. 

First, purely attractive potentials, gV < 0, where the 
two and three-body thresholds coincide, g — > [6T|10|. Sec- 
ond, potentials with negative or zero net volume but not 
necessarily weak. These potentials have at least one bound 
two-body state for all, even infinitcsimally small, strengths 
5,18,21 . Thus, again the two and three-body thresholds 
are the same and Borromean states cannot exist. For the 
same reasons these potentials can not provide Borromean 
binding also in ID. The third example is the delta shell 
potential, 

gV = h 2 /(2 f id 2 )(X + S(r/c - 1) - \-5(r/d - 1)) , (18) 

with c — > d. Consider a potential with three bound parti- 
cles, i.e. A+ < AT, Ref. j3] proves that if three particles are 
bound with a potential gV , then two particles arc bound 
with a potential \gV , thus the potential 

h 2 /(2^d 2 )^\ + 8(r/c- 1) - \\-5{r/d- 1)), (19) 

with c — > d binds two particles. From the discussion in 
connection with eq. © for two particles, we then know 
that A + < A"' = A_, and as A + — > we therefore find 
A+ = A_. Again, Borromean systems do not exist for this 
potential. 

Surprisingly, even potentials with positive net volume 
do not always provide Borromean binding. First it was 
pointed out in Ref. [5] that the Lennard- Jones potential 
have dimer and trimer thresholds for the same coupling 
constant. Moreover numerical search for three-body bound 
states with potentials of the form 

9 V(r) = (h 2 /2 f ib 2 )[-a 1 e- r2 / b2 + a 2 e~ Ar2 ' h2 ], (20) 

with cti > shows absence of Borromean binding [MTU] . 
This has to be compared with 3D behaviour, where those 



A. G. Volosniev et al: Occurrence conditions for two-dimensional Borromean systems 



5 



potentials always have a regime, when three but not two 
particles are bound [4ll24] . Further extensive numerical in- 
vestigation carried out in Ref. |10j yielded an example of 

potential gV(r) = {h 2 /2pb 2 )[2 e - r2 ^ 2b2s >-5.7e- 2r2 / b2 } that 
can provide Borromean binding. Another example with 
Borromean binding was suggested for 3D in [3] but it is 
also useful in ID and 2D 



gV(r) 



f(r) 



<Cxb 
r > C x b 



(21) 



where /(r) is a positive continuous function that vanishes 
at infinity. We see that if the constant C is sufficiently 
large we obtain Borromean binding for | = < g < 
gC r = 2. This potential can be used in any dimension to 
obtain the largest possible Borromean window (g^/g^)- 
Moreover eq. (|2"Tj) suggests that Borromean systems can 
be obtained for any decay of the potential at infinity. For 
example, potentials with 1/r behaviour at infinity, that 
we do not consider above, and deep enough pocket can 
produce Borromcan binding. 

The overall conclusion is that to get Borromcan bind- 
ing we need potentials with positive net volume which 
leads to finite strength, in contrast to (infinitesimally) 
small. Moreover, the examples we provided above show 
that so far we do not know potentials without an outer 
barrier that can produce Borromean binding. It might 
mean that only potentials with an outer barrier are able 
to provide Borromean states in 2D. 



3.2 Three-body conditions 

To get Borromcan binding we need the three-body wave 
function to be localised in the region where all three parti- 
cles are close to each other. However localization increases 
the kinetic energy. This interplay between kinetic and po- 
tential energies defines the possibility for existence of Bor- 
romean states. We illustrate it by using the hyperspherical 
expansion method, which is efficient for weakly bound sys- 
tems [251. 126] . An upper bound for the energy is found by 
use of just the lowest hyperradial potential 



d 
dp 2 



where E is the energy, m the mass of the particles la- 
beled (1,2,3), Vcs. is the hyperradial potential, and the 
hypcrradius is an average length coordinate defined by 
3p 2 = Y,i<k r ik> where r ik = r, - r k . We also have p 2 = 
3 r i2 + 1^12,3 an d r i2,3 = j(r? 3 + r§ 3 ). The eq. is 
well-known for 3D [25], where 3/4 is replaced by 15/4. 

First, we consider weak purely attractive potentials in 
2D that support only one two-body bound state with en- 
ergy E 2 and root mean square (rms) radius 



R 2 



^-H 2 /(3fiE 2 



(23) 



In this case, the effective potential, V e s, supports only 
two three-body bound states [28U26U14U11TT5] with en- 
ergies £3 = 16.52i?2 and E 3 = 1.27E 2 and rms radii 



i?3 = 0.305i?2 and R3 = 2.55R 2 . Those states are in- 
dependent of the details of the interparticle interaction 
for sufficiently small E% and are determined by the long- 
range behaviour of V c g, which is defined by E 2 . It follows 
that these two universal states always exist if A + — > AST. 
Consequently Borromean binding occurs when we have a 
third state with rms radius proportional to the range of 
the interaction. Then two bound three-body states move 
into the continuum (become unbound) with E 2 — > 0, but 
the ground state at smaller distance remains bound. This 
mechanism for Borromean states reflects that the short- 
range part of the potential V e g is responsible. 

Let us now make an explicit division of V e f / into two 
parts V e s(p) = V sh (p) + V\ on (p), where Vj on depends just 
on two-body binding energy and may support the two 
weakly bound three-body states, and V s h(p) depends on 
the details of the interparticle interaction. We illustrate 
this division on Fig. [3] A crude way to estimate the in- 
terplay between kinetic energy and potential energy is by 
using the short-range part V s h in eq. (|2"21 instead of V e g. 
Appearence of the bound state in such an equation is dis- 
cussed in Ref. [20], that is 



dp = 0, 



where the function <pi(0,p) is defined through 



(24) 



p 



^( 1+ r dP 'MQ,P')vM 

2 V Jo VP' 
1 rp 



VP Jo 



dp'p' 3 / 2 M0,p')V sh (pi). (25) 



Eq. [2U directly expresses that a vanishing strength of 
V s h cannot satisfy this equation. Thus, since a Borromean 
system requires a bound state, the potential must have a 
finite strength. It means, for example, that square wells 
from eq. (fT4"f with R s /Ri — > can not provide Borromean 
binding because they need a very small depth of the pocket 
to provide a two-body bound state. 



3.3 Square well potentials 

Eq. (|2~4"]l qualitatively expresses the idea of the interplay 
between kinetic energy and interaction. A deep enough 
attraction inside a repulsive barrier, e.g. eq. (|9|), allows 
Borromcan systems, since all three particles can benefit si- 
multaneously from the attraction, without extending into 
the barrier. This is discussed for a schematic case in the 
Appendix. One would expect that potentials without an 
outside barrier, e.g. eq. (|14p . also must support Borromean 
binding for some \-,A±,R s /Ri. However extensive nu- 
merical searches did not reveal Borromean binding. Unfor- 
tunately, exploiting numerical search, we can not exclude 
it rigorously for all potentials without an outer barrier. 
However, using a rough estimate obtained from eq. (1241) 
we can suggest a region where Borromean system might 
occur. To do so we establish a qualitative connection be- 
tween a two-body square-well potential with infinite core 
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Fig. 3. Division of the lowest effective adiabatic potential 
V c ff (solid line) for £2 — > into V B h (left) and V\ on (right). 
Range and strength of V s h are defined through the range b 
and strength g of the two-body potential. Range and strength 
of Vion are universal and defined through the two-body en- 
ergy £2. V[ n supports two universal three body states (dotted 
and dash-dotted lines). Here V s h supports a third bound state 
(dashed) which becomes Borromean for E2 = 0. If V sn cannot 
support a bound state for E2 = 0, a Borromean state does not 
occur. 
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Fig. 4. Axis and the solid red line is the same as on Fig. [2] 
Dashed and dot-dashed blue lines give an insight into the re- 
gions where three-body bound states with root mean square 
radius proportional to the range of the potential (see the Ap- 
pendix) might exist. The space between solid and dashed (dot- 
dashed) line (when solid line is higher) gives a rough estimate 
for h and s that can support Borromean states. 



and an effective three-body potential. We want to use only 
simple solvable potentials to suggest possibilities that un- 
ravel the general trend (see the appendix). 

Using this connection wo determine parameters of the 
two-body potentials that give three-body bound states 
with rms radius proportional to the range of the poten- 
tial. To determine whether this is a Borromean state, we 
compare parameters of this potential with a two-body po- 
tential, that satisfies eq. (fTT)) . We present the result of 
this qualitative investigation on Fig. 01 where two-body 
bound states exist above the solid line, while dashed and 
dot-dashed lines represent parameters of the two-body po- 
tential that might allow three-body ground state with rms 
radius proportional to the range of the potential (see the 
Appendix). Thus Borromean systems are likely to exist 
below the two-body and above the three-body curves. The 
optimum configuration occurs when the two-body attrac- 
tion is most effective for all pairs. By simple geometric con- 
sidetaions, this dictates that Ri must be larger that 2R S . 
On the other hand R s /Ri can not go to zero, where as fol- 
lows from Fig. S] a two-body system is easily bound. This 
analysis suggests that the region where Ri ~ 2R S has to 
produce the largest Borromean window, because the neg- 
ative part of the potential is deep enough to overcome the 
kinetic energy term, and radii do not deviate too much 
from each other, which allows the most effective attrac- 
tion. This qualitative procedure unravels regions where 
three body Borromean state are most likely to exist. 

3.4 Numerical illustrations 

Our numerical procedure is based on the stochastic vari- 
ational method with basis of correlated Gaussians. This 



procedure was proven to give accurate results |27j . how- 
ever the convergence slows down for potentials with a 
hard core. Numerically convenient potentials are linear 
combinations of gaussians, because it allows calculation of 
matrix elements analytically. Moreover extensive numer- 
ical search with different potentials did not reveal Bor- 
romean states without an outer barrier. Thus, we choose 
a potential with outer barrier in the form V± = — ^ — L , 
with F = 2^(exp(-0.5r 2 /& 2 ) - 2 exp(-2r 2 /6 2 )) , where 
the net volume is zero. This potential is attractive for 
r 2 /b 2 < § In 2 and repulsive otherwise. We show in Fig. [5] 
the critical values for binding the two- and three-body 
systems as function of the attractive strength, A_. The 
straight line for the delta-shell potential is followed for 
small A_ for both two and three particles. The zero net 
volume behavior of A^ ps A_ is found as expected in this 
region. All three curves begin to deviate around A_ « 0.7, 
and the finite range potentials reveal their divergent char- 
acter. We define the critical values, X c f and A C Z , above 
which two and three-body bound states respectively are 
present independent of A + . This limit can be described 
by potentials of the form eg. (|21[) and we see numerically 
the expected ratio, A^/X^ - 2/3. When A_ > Ai r the 
attractive pocket alone supports a two-body bound state. 

At the threshold A+ = AST two three-body states dis- 
appear into the continuum. They have the same threshold 
as the two-body system in complete analogy to weakly 
bound systems where the two and three-body thresholds 
are identical [TUll2"5] . The difference for the deep potentials 
is that the Borromean state remains. 
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Fig. 5. The two- and three-body critical strengths, A+ r (dashed 
line) and A c + (dot-dashed line), on the positive part of the po- 
tential as function of the strength on the negative part. The 
potential is gV = X+V+ + A_V_ where V± are positive and 
negative parts of ( exp(-0.5r 2 /6 2 ) - 2 exp(-2r 2 /6 2 )) , re- 
spectively. Solid line represents results for delta-shell potential, 
i.e. A^_ r = A_. 

4 Experimental realization 

From the discussion above it is obvious that Borromean 
systems in 2D are not as abundant as in 3D. In particular, 
we need an outer barrier and an attractive pocket. Sys- 
tems where there is no outer barrier and those with non- 
positive net volume integral of the interaction are there- 
fore out of the picture. This includes the case of neutral 
atoms in a plane which would have van der Waals-type 
attractive pockets and inner hard-core repulsion. A pop- 
ular way to achieve long-range interactions is to use polar 
molecules [29)130] . In particular, layered systems with po- 
lar molcules have promising properties and have recently 
been stabilized experimentally [31] . The geometric dipolc- 
dipolc interaction of polar molecules in a multilayer sys- 
tem supports a number of interesting few-body states [221 
32, 33 , 34 , 35 , 36 that arc indicators of non-trivial many- 
body pairing [3"?]f3"51[3lJ] . However, intcrlayer interactions 
always have zero net volume integral [22], while the in- 
tralayer interactions are purely repulsive [3 5) and do not 
have the necessary pocket plus barrier structure. 

In the following we list possible experimental setups 
where it could be possible to simultaneously maintain an 
outer barrier and a substantial inner attractive pocket. 

The first option is to use cold ionized atoms in a two- 
dimensional geometry. Ideally this would be ions confined 
to live on a surface. Other interesting setups could be 
ions trapped near a surface |40j or in a Penning trap |41| . 
While these studies are mostly done with atoms in a crys- 
tal state, the few-body dynamics studied in the current 
paper requires that the ions also have motional degrees 
of freedom that are continuous in some range (or at least 
quasi-continuous if there is some weak confinement in the 
2D plane of motion). To reach the Borromean regime, we 
have the outer barrier from the Coulomb repulsion of the 
ions. The inner attractive pocket would then need to be 



provided by short-range van der Waals interactions. This 
pocket is then required to sit at some small yet not too 
small distance in order to avoid the regime where every- 
thing is controlled by chemical reaction dynamics. What 
is particularly complicated about this proposal is the fact 
that the outer repulsive Coulomb barrier could be ex- 
tremely large and render any inner attractive pocket irrel- 
evant. This could possibly be counteracted through elec- 
tron screening that lowers the Coulomb barrier. Here we 
imagine that a plasma of ions and electrons could be use- 
ful if it can be confined to 2D. A more realistic possiblity 
could be mobile impurities in solid-state system where the 
background electron density provides a screening effect. 

A second, and presumably much easier, option is to 
use polar molecules in external fields. It has been pre- 
dicted that applications of AC and DC fields in systems 
with particles that have non-zero permanent electric or 
magnetic dipole moments provides a way to taylor the 
inter-molecular interactions [32,43,44,45]. In a squeezed 
geometry (quasi-2D), AC fields can be used to control 
the existance of two-body bound states for both bosonic 
and fcrmionic molecules, and correspondingly may result 
in bound three-body states of AC field dressed molecules 
[4"tj] . What is needed here is the presence of both outer bar- 
rier and inner attractive pocket. As is discussed in Ref. [43] 
this can be achieved by applying both an AC and a DC 
external field to the system. DC fields have been used to 
align polar molecules and control the overall magnitude 
of the dipole moment in two-dimensional systems |31j and 
very recently the application of AC fields to the same ge- 
ometry has been reported [47] . A combination of these two 
external influences could provide the potential profile nec- 
essary to produce and observe Borromean bound states in 
a 2D setup. 



5 Conclusions 

Two-body binding is easier achieved in 2D compared to 
3D, because of the negative centrifugal barrier in the ra- 
dial two-body Schrodingcr equation. In all other spatial 
dimensions the centrifugal barrier is zero or positive. This 
feature makes 2D geometry a case of special interest. 

We investigate possible few-body structures in 2D with 
the goal to find conditions for Borromean states. Thus, we 
search for the potentials that can support three, but not 
two-body bound states. Our focus is on potentials with 
positive and negative parts, as purely attractive potentials 
in 2D always support two-body bound state thus exclud- 
ing Borromean states. 

We show that the necessary condition for a potential 
to provide Borromean states is a substantial repulsive and 
attractive part. This a consequence of the fact that near a 
two-body threshold a Borromean state is localised at small 
radii. This can happen only for potentials with substantial 
attraction to outweigh the kinetic energy of localisation. 
Moreover, our numerical search indicates that potentials 
without outer barrier are highly unlikely to support three- 
body states without a two-body state. 
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We conclude that experimental observation of Bor- 
romean systems in 2D is possible only for interactions 
that have an outer barrier, a substantial attractive re- 
gion, and positive net volume. Polar molecules or ions 
in squeezed geometries are potential candidate systems 
for the occurence of low-dimensional Borromean states. 
A more speculative possibility is to look for Borromean 
signatures among impurities on a solid-state surface. The 
details of the experimental scenarios and which options 
are more viable goes beyond the current principle discus- 
sion and will be the focus of future studies. 
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A Connections to the Hyperspherical 
Expansion 

Here we qualitatively discuss the connection between an 
intcrparticle interaction and and the lowest adiabatic po- 
tential in the hyperspherical expansion method. 

The hypcrradial potential for three bosons interacting 
pairwisc through a square well of radius Rq and depth 
Vo is determined semi-analytically for s-waves by solving 
trancendcntal algebraic equations [48] . The hyperradius, 
p, is defined through two-body distances between pairs of 
the three particles which implies that the 3-D-procedurc 
is directly applicable for 2D as well. The method divides 
the p-coordinate into 4 intervals, that is 

(i) with p < Rq/V2 : where all three two-particle dis- 
tances are smaller than Rq. 

(ii) Ro/y/2 < p < Rq^/2/3, where at least 2 pairs 
interact, but configurations of 3 interacting pairs are also 
possible. 

(iii) Rq^/ 2/3 < p < Rq\/2, where at most 1 pair inter- 
acts, but configurations of 3 non-interacting pairs arc also 
possible. 

(iv) Rq^/2 < p, where at least 2 pairs don't interact. 
The hyperradial attractive potential in turn increases 

from 3Vq for p = through various continuous steps to 
a value less than Vo at p — RqV2. For larger p- values 
the potential approaches zero as 1/p 2 until the scattering 
length is reached for p x a > JJ . For p larger than a 
the potential approaches zero even faster, that is in 3D as 
oc a/ p 3 and in 2D as oc a/(p 3 ln(p)). These potentials must 
be supplemented by the centrifugal barrier term, which in 
3D and 2D corresponds to effective angular momenta of 
I* =3/2 and I* = 1/2, respectively. 

The large-distance behavior in 3D allows a number 
of bound states provided the scattering length is suffi- 
ciently large. This means that it is possible to vary Ro and 
Vo, where the scattering length is maintained to be much 
larger than Rq while the volume is too small to support 
a bound two-body state. By increasing a towards infin- 
ity, still for an unbound two-body system, the number of 
three-body bound states increase logarithmically towards 
infinity. This is the Efimov effect. 



In 2D any, even infinitesimal, overall attractive poten- 
tial provides a bound state. Thus, when the three-body po- 
tential can bind with the I* = 1/2-term, also the two-body 
system with I* = —1/2-term is bound. No Borromean sys- 
tem can be constructed, and no Efimov effect exists, see 
subsection 13.21 

We consider now the more general two-body potential 
of an attractive square well of radius R s and depth V s , 
and a repulsive barrier between R s and Ri of height V/. 
The most interesting combination in the present context 
is when Ri > 2R S , since then the large-distance configura- 
tions apply for the s-part of the potential before the short- 
distance behavior of the ^-part has ceased to be present. In 
other words, the attraction contributes fully for all pairs 
when p < R s />/2. 

The volume, h 3 = k 2 (Rf — R%), for the three-body 
potential is in general much larger than, h 2 , for the two- 
body potential. We can roughly relate by h 3 ps 
where = m/fi = 2 is the ratio of the mass, m, in ea. (|22l) . 
fv = (V)/Vq ps 2 is the average three-body potential in 
units of Vo, and /r ps (Rq^/2/Rq) 2 = 2. Increase of the 
volume leads to the possibility of the Borromean binding. 

In the limit of a huge barrier, Vi is very large, we can 
then evaluate the condition for a bound three-body state. 
The hyperradial wave function, Ji=x(K,p), for I* = 3/2 
must have a node at p = R s /y/2, where k is the corre- 
sponding wave number. The dimcnsionlcss volume ratio 
is then, h 3 /h 2 = fufvfn = 2x3x1/2 = 3. The node, 
14.66, is smaller than 3 times the corresponding node, 5.78 
of the two-body wave function, Ji=o(kr), which imme- 
dately implies that there is a window where three, but not 
two, particles are bound. Borromean systems exist when 
4.89 < h 2 < 5.78, that is for small s < 0.5. 

Changing the sign of the two-body potential gives a 
repulsive core up to R s and afterwards until Ri an attrac- 
tive well. We use again the relation between two and three- 
body potentials where three-body binding must arise from 
the attractive part in the interval R S V2 <p< R1V2. The 
limit of s — > correspond to two-body potentials resem- 
bling an overall very weakly attractive potential. In this 
case Borromean systems cannot exist, see subsection 13.21 

Assuming an average potential strength of 2Vi we get, 
with S3 = s<i = R s /Ri that Borromean systems exist when 
8/12 > ft.3. This comparison can be improved by replacing 
2Vi by a weighted average of the attractive strength, (V) ps 
JVrdr/ J rdr = 19(l-36s 2 /19)/[12(l-s 2 )]. The average 
decreases with s to a value of 10/9 for s = 1/2. This 
curve, shown in Fig|4] - dot-dashed line, allows a rather 
large window for Borromean states in an interval from 
O.K s < 0.5. 

These estimates of restrictions are rather conservative. 
First they are only estimates, and second two other, so far 
neglected effects, oppose existence of Borromean states 
derived from the hyperspherical formalism. The first is 
that the approximation to use only the lowest adiabatic 
potential implies that the correct three-body energy must 
be lower. This is derived from the general theorem that 
the correct energy lie between the results from the lowest 
adiabatic potential with and without the diagonal non- 
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adiabatic term. This tends to further narrow down the 
Borromean window. 

The second effect is that the fraction of the 2D coordi- 
nate space where the distance between two particles is less 
than R s is R 2 s /(2p 2 ) for a given hypcrradius. Then one of 
the attractions, Vi, should be replaced by the repulsion, V s , 
in total amounting to V] — > VJ(1 — s 2 (l + V s /Vi)/3) when 
p = R1/V2. For a larger p = Riy/2/3, the factor 1/3 
should be replaced by 3/8. The requirement of an over- 
all positive volume of the two-body potential provide the 
inequality s (1 + V s /Vi) > 1, which gives a limit on the re- 
duction factor, that is at least (l-s 2 (l + y s /Vi)/3) < 2/3. 
Thus, we should further restrict the volumes by multiply- 
ing with 2/3. This curve is also shown in FigQ] (dashed 
line), leaving still room for Borromean states. However, 
a reduction by an additional factor of 2 would close that 
window completely. 

In conclusion, the features of the square well must be 
maintained in other potentials, that is with both attrac- 
tive and repulsive parts. The conclusions are therefore 
much more general. Whether a small window is open for 
Borromean states in 2D without a confining outer bar- 
rier remains to be seen. The features of candidate po- 
tentials are that the ratio of attractive to repulsive vol- 
umes must be between 8 and 15. The ratio of repulsive 
to attractive strengths must be limited to be smaller than 
1/s 2 - 1 < V s /Vi < 2/s 2 -\. 
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